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This article describes firstly a wind power production line, principally a wind 
turbine constitutes her and brushless doubly fed induction generator 
(BDFIG). The models of these components are developed, and control 
objective of BDFIG is to achieve a dynamic performance similar to the 
doubly fed induction generator (DFIG) using a stator flux field oriented 
control (FOC) and direct power control (DPC) strategy. After, the simulation 
program using Matlab/Simulink has been developed. The performances of 
this strategy are evaluated and analyzed, so the results shows a good 
robustness great dynamic and a precision of speed. 
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1. INTRODUCTION 

Wind power has become increasingly popular because of the increasing difficulty of the 
environmental pollution and the greenhouse effect. Huge efforts have been made in promoting the wind 
energy conversion system , to reduce costs, increase reliability and robustness [1]. The permanent magnet 
synchronous generator (PMSG) and the doubly fed induction generator (DFIG) have become the most 
popular generators for wind energy applications. On one hand, the use of brushes and slip rings associated 
with the rotor of DFIG decreases the robustness of system and increases the maintenance cost [2]. 

The cost of maintenance for traditional DFIG based windgenerators increased the pressure to seek 
other alternative generator systems. TheBrushless Doubly-Fed Induction Generator (BDFIG), also known as 
a self-cascaded generator, is composed of two stator of different pole numbers called stator of power winding 
(PW) and stator of control winding (CW) and a special rotor winding. Normally, the two stator supplies are 
of various frequencies, one a fixed frequency supply linked to the grid via converter, and the other a variable 
frequency supply derived from (AC-DC-AC) converter. The natural synchronous speed of the machine is 
equal to: 


_ a) p +C °c 

P p +P c 


( 1 ) 


Where P p , P c , (x ) p , 0) c and 0) r represents respectively, pole-pair of the PW, pole pair of CW, angular 
frequency of PW, angular frequency of CW and , Then the shaft angular speed. 
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The variable speed constant frequency (VSCF) generator system of the turbine is more important to 
improve the effectiveness in capturing the maximum energy of the wind and the assistance of the high 
quality, the efficiency and the power controllable, where the major challenge is the independent control of the 
active and reactive power exchanged between the BDFG and the grid. 

A generator of a wind turbine should be fully controllable so that it can operate at a speed of shaft 
which is dependent on wind conditions to obtain the maximum output power. Furthermore, reactive power 
control is important for the production of electricity, and consequently, the functionality of the generator 
controller to regulate the power factor is essential. Because the BDFIG isn't stable over the complete running 
speed range, a controller is required to stabilize the machine at the same time as achieving nice dynamic 
overall performance in controlling the speed and reactive energy. Robustness and fast dynamic reaction are 
critical capabilities of such a controller for wind-electricity programs [3]. 

Some strategies of control have been used up to now in this machine (scalar current control, direct 
torque control, fuzzy power control, sliding mode power control, and the rotor flux oriented control. A new 
vector controller using a dynamic model with a unified reference frame based on the PW flux was 
investigated for the BDFM [4]. Furthermore simplified controller oriented with the PW stator flux with a 
complete mathematical derivation frame has been exhibited, with some experimental results presented on 
both speed and reactive power regulating [5]. 

Classical control strategies for grid side converter (GSC) systems regulate both the active power and 
reactive power flow by controllingthe currentvector orientation with respect to the grid voltage vector. This 
technique is referred to as voltage oriented control (VOC) [5],[6]. The VOC technique decomposes the AC 
line currents into the direct and quadrature components by using the Park transformation. This makes 
possible to regulate the active and reactive powers by controlling the decoupled AC currents, using current 
controllers [7]. 

However, there are a new strategies used to control grid- tied VSI systems in this type of program 
which include as the Direct Power Control (DPC). One of the fundamental qualities of the DPC is that the 
control directly the instantaneous active and reactive powers instead of instantaneous AC line currents. The 
DPC is mainly based on the principle of Direct Torque Control (DTC) for electrical machines [8]. In effect 
the DTC guides the stator flux and controls the torque of AC machines according to a switching table. 
Similarly, the DPC regulates instantaneous active and reactive powers by using an optimal switching table. 
This table determines the VSI switchingstates by means of using the errors of power and the position of the 
grid voltage vector. 

This paper is dedicated to the study and control of the wind power system based on a BDFIG the use 
of indirect control power of machine side converter (MSC) by using PI controller and (DPC) control of 
(GSC) converter to control the voltage of the DC Link, MPPT method is carried out to extract the maximum 
power available. Figure 1 show the schematic diagram of power and control circuits of WECS studied in this 
paper. 



Figure 1. Schematic diagram of a BDFIG based wind power system 
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2. WECS MODELLING 
2.1. Aerodynamic Model 

The wind turbine is a device that converts kinetic energy of wind into mechanical energy. Wind 
power is defined as follows [9]: 

P v =^pR 2 v 3 (2) 

The aerodynamic power ( P a ) captured by the wind turbine is given by: 

P a =^K 9 R 2 C p (k,$y (3) 

The wind turbine is characterized by its curve Cp=f(A) [10]-[13].The tip speed ratio (k) is expressed by the 
following expression : 

R£l t 

X =-— (4) 

v 

C p : power coefficient; p\ density of air (1.25 Kg/m3); 
v: wind speed; W tur : wind turbine speed; R : turbine rayon. 

The power coefficient C p represents the aerodynamic efficiency of the wind turbine, it is determine 

by [10]: 


Cp(l , Z?) — Cj(c 2 // ■ + c 3 b + c^)c 5 + Cg/ ■ 


(5) 


Where, 


l t l + 0.08Z? b 3 + 1 

With, 

c, = 0.51,c 2 = 116,c 3 = 0.5,c 4 = 5,c 5 = 21 and c 6 = 0.0068. 

Figure 2 show the typical curves of the power coefficient. We note that the maximum (C pmax ) is 
reached for: (/ = 8 and b =0). 



Figure 2. Typical curve of power coefficient 
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From this power the wind torque is given by: 


r p a r r - Sy * 1 

C = -= Cp. -.- 

a W 2 W 


The dynamics equation is : 


J.^=C 

dt 


J =-^+/ 


(7) 

( 8 ) 

(9) 


With: 

C 

mec : Total mechanical torque applied to the rotor; 

J : is the total inertia that appears on the generator rotor. 
The mechanical torque: 


c = c - c - c 

mec g em vis 


( 10 ) 


2.2. BDFG model 

The model of the brushless doubly fed induction generator in the PW synchronously rotating d-q 
reference frame is expressed as [14]-[16]. 


JO„ 

y _ D • _|_ 
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V dsp ~ R sp^dsp + " 


JO 


dsp 


dt 


^d sc R sddsc ^ 


JO 


dsc 


oo n O„ 




qsc 


Vqsc R sdqsc ^ 


dt 
JO 

—+K - (Pp + Pc K )®dsc 


0 - R ddr + 


JO 


dr 


dt 
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^ dsc ^sddsc " 1 " scddr 

F =L i + L i 

qsc sc qsc scr qr 

F dr = L r idr + ^ scddsc P spddsp 

F —Li + AL i + L i 

qr r qr scr qsc spr qsp 


( 12 ) 


The electromagnetic torque is expressed as [10] : 

C e m — PpM S pr(iqspidr~idspiqr )~^Pc^scr(^dsc^qr~^qscLr ) 

The active and reactive powers of the stator of PW are defined as: 


(13) 
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^P 2 K 'Isdp Vgqp -I sqp ) 


(14) 



(15) 


3. CONTROL STRATEGIES 

3.1. Maximum Power Point Tracking Technique 

It is evident that it is clear that the wind power is affected by the wind speed. The wind speed 
increases with the height most rapidly near the ground, increasing less rapidly with greater height.The wind 
speed at which electric power production starts called the cut-in wind speed. The turbine will developenough 
mechanical power to rotate itself at slightly lower speeds, but this wind speed will actually supply all the 
generator and transmission losses so that useful electric power cannot be produced. At rated wind speed the 
power input to the wind turbine will reach the limit for continuous operation (rated power). When the wind 
speed exceeds this level the excess power in the wind must be discarded by varying the pitch angle of the 
blades to prevent the turbine overloading. The power is maintained at its rated value until a maximum wind 
speed is reached the cut-off wind speed (Vcut-off) then the turbine will shut down Figure 3 [13]. 



Figure 3. Wind turbine power curve characteristics 


Maximum power point tracking (MPPT) strategies play an important role in wind power conversion 
systems (WECS) because they maximize the power extracted from thewind, and therefore optimize the 
conversion efficiency [17]. 

Figure 4 shows the Power-Speed characteristics of the wind turbine, the peak power for each wind 
speed occurs at the point where Cp is maximized. To maximize the power generated, it is therefore desirable 
for the generator to have a power characteristic that will follow the maximum C pmax line. 
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Figure 4. Characteristics curve of wind turbine 


To extract the maximum power generated, we must maintain X at the optimal command rotor speed 
X opt . The coordinates of the optimal point is the maximum power coefficient Cp are ( A, opt =8.1, C pmax = 0.48, 
P = 0). Increasing p allows the reduction of mechanical power recovered from the axis of the wind turbine 
(see Figure 4). 

Two strategies are used in literature [18], with or without speed control. In this paper, we used the 
strategy with speed control; it permits to conduct the speed wind turbine to the desired value which 
corresponds to the maximum power point. The simplified representation of wind turbine model with speed 
control in the form of diagram blocks is given in Figure 5. 




Figure. 5 Wind turbine model with speed control 


3.2. Control of the BDFG with a PW field oriented 

If the d-axis of the PW synchronous reference frame is aligned with the PW air gap flux, the R p of 
PW is neglected.Then, the relation between the PW voltage and its flux is: 


f^ = o 

\K,p = V P = w p f p 


(16) 
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| ^ dsp L sp ifap "P A/ S p r i^ r 

1 0 = L i + M i 

J sp qsp spr qr 

From (17), the equations linking the rotor currents to the PW currents are deduced as: 


l 



M ct 


(17) 


(18) 


Figure 6 shows the control algorithm of BDFIG 



Figure 6. General Scheme control of BDFIG 


This scheme is based on the cascade regulation method [4]. Two independent regulation paths are 
implemented: 


Reactive power control: Q 


V ® La ® V* 


Active power control: P sp ® i spq ® i scq ® V scq 


3.3. PW current Control 

The current derivative (CW) is given by the equation: 


dL 


dt M spr M scr 


R sphp . R s P L r s sp i scd 

' scd 


R 


M w M S cr dt M scr M scr ' 


R sp sp • .q 

——— z + W„l q 


(19) 


di c . 


R spL sp 


R sp sp di SC q 


R, 


dt M spr M scr M spr M scr dt M spr M scr 


sp I" w Rsp ' 


R sp L r s sp 




l spq + W Rsp l scq 


( 20 ) 


Relationships (19) and (20) can be rearranged into two terms: 


di, 


scd 


dt 

di 


scq 


dt 


a xd -hpd a yd (*' spq ’ hcq ’ f sp |) 
—a xq m hpq " 1 " ^yq (hpd ’ hcd ’ sp |) 


( 21 ) 
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With: 
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L Ls 
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The a xd ( i d ) and a xq ( / ) reflect a linear relationship between direct current vector components of 

(PW) and (CW), while ayd( i spq , i scq , |F sp | ) and ayq( i spq , i scd , i scq , |F sp | ) represent the current coupling 
between the cross vector components d and q. 

3.4. CW current control 

The current derivative (PW) is given by the equation: 


di sp d _ M S p } M SC r . di scd L sp R r R 


spd iri spr 1Y1 scr ,,, ^ scd ^sp iV r . , c , , 1V1 spr 1Y1 scr , T . . > 

L7~~t ;- ( M scr- — - -^ l spd+—?sp )+— --KA S sp l spq- W Rsp M scrLcq ) 


dt L Ls 
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dt M spr * M spr 


di Spq _ M spr M scr n , d L q kp R r. R, I. I, M spr M spr , , .. . . 

dt L Ls xr ' dt~ M l ” M L Ls •('* V* + "V 

sp r sp spr spr sp r sp 


The tension of (CW) given by the equation: 


| Vs cd ^xd spd ^yd ( ^spd ’ ^ spq ’ f scq ’ p sp | ) 

1 II 

1 v scq =b xq .i spq + bj i spd ,i spq ,id,\F sp \) 


whose individual terms are: 


l h —r i A-tr l/ • d" 

L r s s P dt 

Y.fr _ M spr M scr R r . M scr R r I I M scr L sp . ^scr / )R S/ 

:: * L Ls spd L Ls ^ Pp r M spq <«Ls Rsp Rsp Pc r sc V cq 

■ ] sp r sp sp r sp spr g^r sp ^ 

:: m 2 di 

\\b =R i +(L - )-^- 

:: * 5C<? Ls sp dt 

L _M w M xr R ri M scr ^ sp - 1 __ 1 k .... M ^L SP; t 

yq T T l spq~ Ayf § ’Pp^RspW sp | + Pp W r ,, l spq + W Rsp ~ Rsp ~ Pc W r) sc\1 l sc 

" L s P V sp M spr§ S sp q M scr ^r S sp Q 

3.5. Active and reactive power Control 

The expression of reactive power based on the flow and current is: 


Q =-.w F i ,- 

^sp 2 s p s p s P d dt spq 
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3 i 
P = — .w 

sp 2 s p I 


3 d F s 


2 dt 


-■hpd + 2r 


(30) 


A s ( — sp =0 ) the term can be neglected in (i spq ). The transfer function of the control Q sp must consider the 

dynamics of the closed loop control of ( i spd ): 


hjo 


= ^- W spK 


■FT(s l 


The transfer functions for regulating P sp is: 


rjs) 3 , 

— 5 --= —.W F 

hjs) 2 Ip| 


■FT(s) R , gi 


(31) 


(32) 


3.6. Direct power control of GSC 

The basic principle of the DPC is based on the well know Direct Torque Control (DTC) for 
induction machines. In the DPC, the active and reactive powers replace the torque and flux amplitude used as 
the controlled output in the DTC [8],[19],[20]. Figure 7 shows the configuration of the direct instantaneous 
active and reactive power controller. Where the hysteresis comparators and a switching table are used. 



Comparing the instantaneous active power (P) at a reference power (P ref ), this latter is obtained by 
the DC voltage control block at the capacitor terminals, where we use a PI controller (Proportional, 
integrator) tocontrol the error between the sensed voltage (continuous) and referenced8],[20],[21]. Whilst to 
achieve a unity power factor, reactive power reference is directly imposed zero (Q re f =0). 
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Figure 8. Twelve (12) sectors on stationary coordinates to specify voltage vector phase 


For this purpose, the stationary coordinates are divided into 12 sectors, as shown in Figure 8, and the 
sectors can be numerically expressed as: 

(n- 2)* — £ q n £ (n- 1 )*-;n= 1,2,3.... 12 (33) 

6 6 


= |1 forPi P ref - H P 
P 1 0 forP 3 F re/ - H P 


(34) 


Similarly for reactive power controller: 

d JUorQ£ Q ref - H q 
Q | OforP 3 Q ref - H q 


(35) 


With, d v ,d q are digitized error Signals of P and Q ; 0n: voltage phase ;(S a , S b , and S c ) : switching state. 


Table 1. Switching Table for DPC 


dp 

d Q 

01 

02 

03 

04 

05 

06 

01 

08 

09 

010 

011 

012 

1 

0 

101 

111 

100 

000 

110 

111 

010 

000 

Oil 

111 

001 

000 

1 

1 

111 

111 

000 

000 

111 

111 

000 

000 

111 

111 

000 

000 

0 

0 

101 

100 

100 

110 

110 

010 

010 

Oil 

Oil 

001 

001 

101 

0 

1 

100 

110 

110 

010 

010 

Oil 

Oil 

001 

001 

101 

101 

100 


The instantaneous input active and reactive powers of three phase rectifier are generally defined as: 


e a-i a + e b-h + eJc 

(36) 

i e b - e c X + ( e c - e a X + ( e a - e „ X jj 

(37) 



arctang—± 

(38) 

& 
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4. SIMULATION RESULTS 

The simulation was performed with the Matlab / Simulink software to validate the commands 
studied in this work.We realized two operating points, one on the step of speed (see Figure 9(a))and other 
variable speed (see wind profil Figure 16(a)).The wind speed is modelled by deterministic form a sum of 
several harmonics: 


Kind= v o+ a v iSin(Wjt) (39) 

i=0 

The parameters of BDFIG used of this study are given in Table 1. 


Table 3. BDFIG Parametres [4] 


Power Windin (PW) 

Control Windin(CW) 


Rotor 

R,„ = 1-732 w 

R sc = 1 079 W 

R r = 

0.473 w 

L sp ^ 714.8m// 

L c = 121.7m// 

L = 

132.6 mH 

M spr = 242.1m// 

M spr = 59.8m// 

f P= 

:50HZ 

P„= 1 

P c = 3 



=142mH 

C =1.65 O 




Concerning the results of the simulations, two cases have been analyzed corresponding to the 
profiles of wind respectively represented , by the Figures 9(a) and 16(a). However for the two (2) profiles : 

The speeds (see Figures 9(b) and 16(b)) are almost adapted to the wind speed, respective, resulting 
in a very significant increase of the power. By elsewhere, the speed is less than that of the synchrony (750 
tr/m). This technique of extraction of the maximum of power is to determine the speed of the turbine which 
allows you to get the maximum of power generated. Also, we note that the amplitudes of the current of 
phases as the winding of power and control (Figures 10(a), 11(a), 17(a) and 18(a)), consequencement change, 
with the variation of the speed of the wind and that their forms sinusoidal are given by the figures (10(b), 
11(b), 17(b) and 18(b)). 

First According to the Figures( 12(a), 19(a)) and their zoom illustrat by 12(b) and 19(b), 
respectively, it may be noted the robustness of the vector command in term of decoupling and the good 
results obtained by the regulation the active and reactive power, second, we show that a small variation of 
the wind can induce a large variation in the power extracted (mechanical power), because of the 
proportionality of this last to the average value of the cubic speed of the wind. 

We distinguishes the sinusoidal shape of the current of phase of the PW and the reactive power is 
zero, which translated the production of electrical energy under unit power factor. 

The MPPT strategy we allows to provide the totality of the active power produced at the electrical 
network with a unit power factor. However, the command decoupled from the active and reactive power 
allows to regulate the active and reactive power provided to the network. Figures 13 and 20 shows the CW 
current components idc and iqc. This components imposent, respectivement, the reactive power (Q) and 
active power (P). 

The component in squaring of the flow (fluxqp) of (PW) cancels in permanent regime, and the direct 
component of the flow (PW) (fluxdp) is equal to the value (-1.1 Wb) (see Figure 14 and 21), and the brusque 
variation of the active power (Psp) has no influence on the flow of (PW). The Figures 14 shows the flow of 
BP, one distinguishes the orientation of the flow on the direct axis, such as fluxdp = 0, this justifies that the 
command to flow oriented. 

The variation of the electrical power delivered to the network is adapted to the variation 
consequence of the speed of the BDFIG, and the latter is adapted to the variation in the wind speed of the 
wind. This shows the influence of the variation of the mechanical speed as a function of the speed of the 
wind on the electrical power produced. 

It was also noted that the simulation results show a continued good set point for the stator currents 
(PW) and (CW). The instructions of power are well followed by the generator as well for the active power 
that for the reactive power which is maintained zero. The stator frequency of the machine (control winding) 
depends on the speed of rotation of the turbine. 
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Figures 15 and 20 shows that the voltage of continuous bus (Vdc) stabilizes on its reference value 
(600 V) imposed by the command, and The continuing tension in output of the rectifier is well controlled 
and quasi-insensitive to variations in the speed. 
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Figure 9. Wind Speed, Rotor Speed 
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Figure 10. Phase PW current of phase-a and zoom 
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Figure 11. Phase CW current of phase -a and zoom 
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Figure 12. Power of PW 



Figure 13. d-q CW Current 



Figure 14. d-q PW Flux 
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Figure 15. DC link Voltage 




Figure 16. Wind Speed, Rotor Speed 



0 5 10 5.92 5.94 5.96 5.98 6 6 02 6 04 

lime (*) 

Figure 17. Phase PW current 
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Figure 18. Phase CW current 



Figure 19. Power of PW 
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Figure 20. d-q CW Currents 
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Time (i) 

Figure 21. d-q PW Flux 
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Figure 22. DC link Voltage 


5. CONCLUSION 

This work dealt with the modelling and control of a wind system with variable speed based on a 
brushless doubly fed induction generator. First, we have interested in modeling various parts of the wind 
turbine system, in order to establish the various controls on the two converters (MSC and GSC). The PW flux 
vector control is oriented to apply independent control of active power and reactive. Direct power control 
applied to the GSC converter to maintain constant DC link voltage and the technique of MPPT is applied to 
provide all of the active power generated to the grid with unity power factor .The results showed that the 
active and reactive power of the wind energy system could be controlled independently while ensuring 
optimum active power supplied to the grid, DPC applied to the GSC provides excellent dynamic response 
and good steady state performances. The results of the simulation show that performance system is 
encouraging. 
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